Dynamically assisted Schwinger mechanism 
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We study electron-positron pair creation from the Dirac vacuum induced by a strong and slowly 
varying electric field (Schwinger effect) which is superimposed by a weak and rapidly changing 
electromagnetic field (dynamical pair creation). In the sub-critical regime where both mechanisms 
separately are strongly suppressed, their combined impact yields a pair creation rate which is dra- 
matically enhanced. Intuitively speaking, the strong electric field lowers the threshold for dynamical 
particle creation - or, alternatively, the fast electromagnetic field generates additional seeds for the 
Schwinger mechanism. These findings could be relevant for planned ultra-high intensity lasers. 

PACS numbers: 12.20.Ds, ll.15.Tk, 11.27.+d. 
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As first realized by Dirac jl[, a consistent relativis- 
tic quantum description of electrons necessarily involves 
negative energy levels, which - in the Dirac-sea picture - 
are filled up in the vacuum state. This entails the strik- 
ing possibility of pulling an electron out of the vacuum 
by means of some external influence, such as a (classi- 
cal) electromagnetic field [2|, where the remaining hole 
in the Dirac sea is then associated with a positron. Of 
course, to create such an electron-positron pair out of 
the vacuum, one has to overcome the energy gap of 2mc'^ 
between the filled and the empty levels. There are ba- 
sically two main mechanisms for doing so: In a strong 
electric field E over a sufficiently long distance L, "vir- 
tual" electron-positron pair fluctuations may gain this 
energy when qEL > 2mc^. This pair creation process 
is called the Schwinger mechanism [3|, |j] and can be 
understood as tunneling through the classically forbid- 
den region (energy gap). Thus it is suppressed exponen- 
tially 0{exp{—TTEs/E}) for weak fields E, where Es — 
rn^c^ /{hq) is the Schwinger critical field. For E ~ E^, 
the work done by separating the electron-positron pair 
over a Compton wavelength is of the order of the energy 
gap 2mc^. Alternatively, a classical time-dependent elec- 
tromagnetic field will also create electron-positron pairs 
in general (dynamical pair creation). However, if the 
frequency uj of the external field is not large enough, 
huj < 2mc^, these non-adiabatic corrections correspond 
to higher-order (i.e., multi-photon) processes and are also 
suppressed exponentially exp{— ©(l/w)} for small lo |5|. 
These pair-production processes are fundamental predic- 
tions of quantum electrodynamics (QED), but only the 
multi-photon production process has so far been observed 
experimentally: the positron data taken at the SLAG E- 
144 experiment have convincingly been explained by n- 
photon production with n ~ 5 [6|. However, a verification 
of the Schwinger mechanism has still remained an exper- 
imental challenge [7|. Since the Schwinger mechanism 
is non-perturbative in the field, its discovery would help 



exploring the non-perturbative realm of quantum field 
theory in a controlled fashion. Here, we propose a new 
mechanism which can help to overcome the strong expo- 
nential suppression. The basic idea is similar in spirit to 
ideas in the study of atomic physics in strong fields, where 
new experimental and theoretical results show that con- 
trolled engineering of special electric field pulse shapes 
can enhance certain interesting physical processes, such 
as high-harmonic generation and above threshold ioniza- 
tion (for reviews see [8|). 

Many previous theoretical studies of pair production 



[ii,[i3,[ii 



13 have been motivated by the seminal work 
of Keldysh [13| on atomic ionization in time-dependent 
electric fields; in particular the crossover between the two 
main mechanisms of pair creation due to strong constant 
electric fields and due to those with spatial or temporal 
variations has been of interest. It turns out that spatial 
variations tend to diminish the pair creation rate [lOL Il2l | 
whereas a time-dependence typically increases the effect 
[111 . Il4{ . However, a realistic experimental situation is 
usually far more complex and may involve various fre- 
quency and amplitude scales over a wide range. This mo- 
tivates us to study electron-positron pair creation in the 
presence of a strong and slow electric field plus weak and 
fast electromagnetic wiggles. We assume that the slow 
electric field E is strong but still far below the Schwinger 
limit Eg , and that the frequency of the weak electromag- 
netic wiggles is smaller than twice the electron mass. As 
explained before, the pair creation rate of each effect sep- 
arately is strongly suppressed in this case. As we shall 
demonstrate below, however, their combined impact may 
be much stronger, i.e., yield an enhanced pair creation 
rate. These findings could be experimentally relevant in 
view of the next-generation light sources 15] aiming at 
approaching the Schwinger limit via high-harmonic fo- 
cusing, which typically generates a high-frequency tail. 

For a first estimate, let us treat the strong and slow 
electric field adiabatically and non-perturbatively (as our 



Hq problem) and the weak and fast electromagnetic wig- 
gles non-adiabatically and perturbatively. In terms of 
the field operator ^ and the Dirac matrices a and /3, the 
unperturbed Hamiltonian density reads {h = c= 1) 



Ho = i''' {ia- V + m(l+ V) "if , 



(1) 



where V{r) = qAq encodes the (approximately static) 
electric field. The perturbation Hamiltonian 



Hi = qi'^a- A{t,r)i' 



(2) 



contains the vector potential A{t, r) of the fast electro- 
magnetic wiggles (scale separation) . The Hq problem can 
be diagonalized via the usual normal-mode expansion 



*(i,r) 



aiui{r)e-"^'* + S|w/(r)e+*"^* 



(3) 



where the positive/negative energy spinor solutions ui{r) 
and vi(r) depend on spin a and wavcnumbcr k, which 
are combined into the index / = {cr, fc}. The elec- 
tron/positron creation/annihilation operators d/, d^ and 
hi,b*j are time- independent. After insertion of ([3]) into 
(PI) , the pair-creation amplitude 21/ j can be calculated in 
first-order perturbation theory 



^ij — 1 / d'^xu]{r)a.- A{t,r)vj{r)e 



iuJit-^-iu jt 



(4) 



Within the present estimate, we are primarily interested 
in the size of the exponent governing the exponential sup- 
pression of the above amplitude. Therefore, we study a 
1-|-1 dimensional toy model, which should reproduce the 
exponent correctly, but disregard the sub-leading pre- 
factor. Let us consider a constant electric field E over 
an interval of length L with vanishing field outside, see 
Fig. [TJ For simplicity, we assume qEL > 2m, i.e., we 
are just above the threshold for the Schwinger effect. In 
terms of the Klein paradox [1^ language, the case lies 
near the border between the intermediate and strong po- 
tential regime. 

For modes /, J whose frequency sum loi + ujj cor- 
responds to the typical frequency of the perturbation 
A(t, r), the exponential suppression of the pair-creation 
amplitude ^ is basically determined by the spatial over- 
lap of the modes u\{r) and vj{r). Assuming uij = toj = 
uj for simplicity (other distributions only induce a shift 
in X but lead to the same result as long as loj +ujj = 2w), 
the spinor ui{r) describes electron modes which are inci- 
dent from the right and totally reflected (due to a; < m) 
by the strong field E, whereas vj{r) corresponds to 
positron modes which are incident from the left and also 
totally reflected. The classical turning points are given 
by x± = ±(m — uj)/{qE), see Fig. [TJ As the electric field 
E is assumed to lie far below the Schwinger limit, we may 
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FIG. 1: Sketch (not to scale) of the level structure (top) and 
the mode functions u] and vj (bottom). The upper and 
lower surface of the Dirac sea at ±m are denoted by solid 
lines, which are distorted by the electric field E in the inter- 
val -L/2 < X < +L/2 with qEL = 2m (top). The horizontal 
dotted lines at ±tj represent the electron/positron levels u] 
and vj with the classical turning points at x±. 



employ the WKB approximation and estimate the expo- 
nential suppression by the integral of the eikonal between 
the classical turning points 



2 

dx ^Jra^ - {qEx - w)2 == ^^ x / ( — ) , (5) 



with fix) = TT -I- 2 arcsin(I - 2x) + 4.^(1 - 2x) ^T^X, 
which can be approximated by fix) ~ 27r(l — x) in the 
relevant interval x g [0, 1]. In the limit x = 0, we exactly 
recover Schwinger's exponential factor exp{—TTEs/E} for 
the pair-creation probability in a static field. For x > 0, 
however, the exponent is reduced. Intuitively speaking, 
the particles do not have to tunnel all the way from —L/2 
to +L/2 because the frequency oj helps them to penetrate 
the strong field region up to x±. In a dual picture, the 
particles tunnel through part of the gap, until they can be 
excited by frequency to. At x = 1/2, we get exactly half 
the exponent, and hence the pair-creation rate is dras- 
tically enhanced exp{— 7rii^s/(2i?)}. Note that, without 
the strong electric background field E, a single photon 
with CO — m/2 would not have enough energy to cre- 
ate an electron-positron pair. Below threshold u < m, 
pair production requires multi-photon processes which 
occur at higher orders in the above expansion scheme. 
For X = Ij the exponent in Eq. (O vanishes as expected, 
since the electromagnetic wiggles have enough energy for 
pair creation to = m. 

The above approach based on the scenario in Fig. [T] 
has the advantage of allowing arbitrary wiggles A(t,r), 
but has the drawback that multi-photon processes re- 
quire high-order calculations. Also realistic backgrounds 
Eit, x) are difficult to handle even though the tunnel- 
ing exponent is expected to be universal for slowly vary- 
ing fields. Inhomogeneous backgrounds as well as multi- 
photon physics can be dealt with by the worldline in- 
stanton technique 1J| which we apply to the following 
specific example: we consider a strong and slow electric 



field pulse superimposed by a weak and fast pulse, botti 
spatially homogeneous, 

E 



E{t) = 



-e,+ 



(6) 



cosh (fit) cosh {ujt) 

with £;s > S > e > and m > w > fi > 0. With 
only one such pulse (say e = 0), the corresponding pair 
creation rate can be computed analytically [10]. For the 
superimposed dual-pulse form in ([5]), we can compute 
the pair creation rate semi-classically using an analytic 
continuation to Euclidean time 0:4 

E £ 

^3(2^4) = ^*7T tan(f7a;4) — i— tan(ti;a;4) , (7) 

\l UJ 

with the tunneling exponent being related to the ac- 
tion of the worldline instanton. Starting with the world- 
line representation of the path integral, we may use the 
electron mass m as a large parameter for the saddle- 
point approximation. The saddle points corresponding 
to the tunneling events are word-line instantons x^ = 
[0,0, a;3(A),X4(A)] which are closed loops in Euclidean 
space-time satisfying the equation of motion 

2 /-Etan(r2a;4) etan(ajX4)\ 



(ia;4 
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1, 



(8) 



where dX^ = dx^ + dx\ is the proper time. This equation 
describes the classical motion of a particle in a potential. 
For small £, the second term tan(ci;a;4) acts as an infinitely 
high rectangular well potential and just reflects instanton 
trajectories Xi{X) at the walls 07x4 — ±7r/2. Between the 
walls, we have an approximately harmonic oscillation due 
to f2 ^ tj and thus i&n{U,Xi) ~ flxi. The structure 
of the solution X4(A) depends on the combined Keldysh 
adiabaticity parameter 

Note that the relevant Keldysh parameter in this multi- 
scale problem is formed out of the dominant frequency lu 
of the fast pulse on the one hand, and the dominant field 
strength E of the slow pulse on the other hand. For small 
7 <C 1, we approach the pure Schwingcr limit, whereas 
large 7 do not correspond to a pure multi-photon regime 
[5| as measured in the SLAC E-144 experiment [a]; large 
7 still involve both multi-photons of frequency to as well 
as a non-perturbative dependence on E. 

For small Keldysh parameters 7 < 7r/2, the instanton 
trajectories do not reach the walls and reflection does 
not occur, i.e., the tan(a;x4)-term has no impact. In this 
case, the weak pulse is too slow to create pairs signifi- 
cantly and we essentially reproduce Schwinger's result. 
Beyond this threshold, 7 > 7r/2, however, the instanton 
trajectories X4(A) change due to reflection at the walls, 
and the instanton action becomes modified 



Anst = m (p dX 



dxi 

Hx 



= Am 



dxi 
dX 



^ (10) 



due to the reflection points at 



rn / TT 

A, = — — arcsm 



qE \2j^ 

Consequently, we obtain (for 7 > tt/2) 



./Vi] 



qE 



2 arcsin ( ^ j + -^^f^f^^ 



(11) 
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FIG. 2: Plots of the instanton action [in units of rri^ /{qE)] 
for the electric field in (|6]), computed using the wordline in- 
stanton method, and plotted as a function of the combined 
Keldysh parameter 7 defined in (|9]). The upper [red] dots 
correspond to uj = 100f2 and E = lOOe, while the lower [blue] 
dots correspond to w = lOfi and E = lOe. The solid lines 
show the Schwinger value of tt, estimated in the text to be 
valid for 7 < tt/2, and the expression (|12p . estimated in the 
text to be valid for 7 > ■k/2. The numerical results agree very 
well with these estimates in the relevant limit where E ^ e 
and oj ;:^ f2. 

At the threshold, 7 = 7r/2, we reproduce the Schwinger 
value -4inst = '^'m-'^/iQE), as one would expect. Above 
the threshold, 7 > 7r/2, the instanton action .Ainst de- 
creases significantly. For example, for 7 = tt, it is re- 
duced by about 40%. For 7 ^ 00, it decays as I/7 
in agreement with the expected multi-photon behav- 
ior [5|. For larger e, the threshold behavior becomes 
smoother, see Fig. [2l Since the pair creation probabil- 
ity, determined by the imaginary part of the effective 
action r[yl^] = — iln (injout), depends exponentially on 
the instanton action, i.e., the saddle-point value 



3(r[A^])^exp{~Anst}: 



(13) 



such a reduction of Ainst implies a drastic enhancement of 
the pair creation probability Q(r[A^]); e.g., a reduction 
of 40% in the exponent could make the difference between 
a suppression of 10""'^° and 10^®, which could mean a 
few electron-positron pairs per day, instead of one per 
year. Of course, one could also reduce -4inst by a factor 
of two via doubling the field E. However, such strong 
fields are at the edge of present experimental capabilities 
and focusing two ultra-high intensity pulses into the same 
space-time region is much harder than superimposing the 



strong pulse and a weak high-frequency pulse. Similarly, 
increasing the characteristic frequency of the ultra-high 
intensity pulse is much harder than just adding a weaker 
field of high frequency. In fact, the envisioned generation 
mechanism (high-harmonic focusing) for the ultra-high 
intensity pulses typically induces a high-frequency tail 
automatically (see below). 

Due to the rather general nature of the arguments used 
above, they are not restricted to the specific pulse in 
Eq. ([6]). One would obtain analogous results for 

A3(X4) ^ -^Ex,f{nx,) - le ^ f , (14) 

with arbitrary positive integers l,m,n and suitable lin- 
ear combinations of such fields. The electric field E{t) 
is then an (even) polynomial in time over the numera- 
tor [1 -|- (wi)^"]'; which might capture some more fea- 
tures of a laser pulse than Eq. ([6]). In this case, the 
threshold is at 7 = 1, but after rescaling ui, we again 
get Eq. (fT2|) . This behavior is confirmed by a world- 
line instanton computation. In order to understand the 
general result, it might be helpful to consider the fol- 
lowing picture: The strong and slow pulse deforms the 
fermionic levels almost adiabatically. As a result, the 
expectation value of the free electron and positron num- 
ber operator (valid for A^ = 0) would scale with E"^ /m^ 
(plus adiabatic corrections E/m? etc.) and thus give 
a rather large result. However, this large number does 
not count real electrons/positrons (e.g., they do not an- 
nihilate) but mostly the instantaneous deformation of 
the ground state. The number of real electron-positron 
pairs, left over after the pulse, scales exponentially with 
exp{— 0(1/_E)}, i.e., non-perturbatively in E [17|. On 
the other hand, the temporary deformation of levels dur- 
ing the strong pulse can be exploited by the small wiggles, 
which can turn the deformed "virtual" pairs into into real 
electron-positron pairs. 

The dramatic enhancement of the exponentially small 
pair-creation probability should be relevant for the 
present experimental efforts aimed at the generation of 
light sources approaching the Schwinger limit, see, e.g., 
[15|. One of the main envisioned mechanisms for the 
final amplification stage is coherent high harmonic fo- 
cusing: Let us imagine sending a laser pulse of ultra- 
high intensity onto a curved metal surface. For a very 
high intensity, the Keldysh parameter of the laser is very 
small and hence electrons in the metal start to oscillate 
coherently and ultra-relativistically. Thereby, they effec- 
tively form a relativistically flying mirror, which reflects 
the incident light with a large Doppler shift and thereby 
generates high harmonics up to large order n. According 
to jlSJ I , the spectrum of these high harmonics is universal 
and the intensity of the n-th harmonic scales with n~^'^ 
up to some cut-off, which is proportional to the third 
power of the Lorentz boost factor of the mirror and thus 
depends on the incident laser intensity. 



For example, an incident optical laser intensity of or- 
der 10^^ W/cm^ reachable in the near future would cor- 
respond to n = 0(10^), which may range up to a signifi- 
cant fraction of the electron mass. Finally, the curvature 
of the metal surface allows us to focus the high harmon- 
ics into a small space-time region, such that the spatially 
and temporally compressed intensity might approach the 
Schwinger limit. In this scenario, the highest harmonics 
(near the cut-off) will not contribute to the peak inten- 
sity significantly - but they still can enhance the pair 
creation probability drastically (compared to what one 
would expect from the Schwinger mechanism alone). 
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